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We present ^®As nuclear magnetic resonance spin-lattice and spin-spin relaxation rate data in 
Ba(Fei_a;Coa;) 2 As 2 and Ba(Fei_a;Cua:) 2 As 2 as a function of temperature, doping and magnetic 
field. The relaxation curves exhibit a broad distribution of relaxation rates, consistent with in¬ 
homogeneous glassy behavior up to 100 K. The doping and temperature response of the width of 
the dynamical heterogeneity is similar to that of the nematic susceptibility measured by elastore- 
sistance measurements. We argue that quenched random fields which couple to the nematic order 
give rise to a nematic glass that is reflected in the spin dynamics. 

PACS numbers: 75.40.Gb, 75.50.Bb, 75.50.Lk, 76.60.-k, 76.60.Es 


I. INTRODUCTION 

The iron arsenide superconductors exhibit multi¬ 
ple phase transitions upon doping, including antifer¬ 
romagnetism, unconventional superconductivity, and 
electronically-driven nematic ordering that breaks Ca ro¬ 
tation symmetry.^ In the context of crystalline materials, 
nematic order refers to an orthorhombic lattice distor¬ 
tion that is driven by electronic rather than structural 
degrees of freedom.^ In the iron pnictides, the transport 
anisotropy far exceeds the orthorhombicity, suggesting 
that the origin is electronic.^ The orthorhombic, or ne¬ 
matic, phase is characterized by the presence of perpen¬ 
dicular twin domains.^ Importantly, there is a strong cou¬ 
pling between the spin and orbital degrees of freedom 
ensuring that the antiferromagnetically ordered Fe spins 
lie along either of these two orthogonal directions.^ Upon 
doping, the nematic and antiferromagnetic ordering tem¬ 
peratures are suppressed, yet strong antiferromagnetic 
fluctuations persist in the paramagnetic state beyond op¬ 
timal doping, even in the absence of long range order.® 
Direct transport measurements of the electronic nematic- 
ity versus strain have uncovered a divergent nematic sus¬ 
ceptibility in the paramagnetic phase. ^ The large nematic 
susceptibility necessarily implies the presence of nematic 
fluctuations in the disordered paramagnetic phase. 

Nuclear magnetic resonance (NMR) has played a cen¬ 
tral role in the investigation of spin fluctuations in the 
iron arsenide superconductors. The ^®As nuclei (/ = 3 / 2 , 
100 % abundant) experience a strong hyperfine coupling 
to the neighboring Fe spins,® thus the spin lattice re¬ 
laxation rate, , is a sensitive probe of the dynam¬ 
ical spin susceptibility of the Fe spins.® In the param¬ 
agnetic state of a homogeneous material, critical spin 
fluctuations exhibit a characteristic time scale, Tc, that 
diverges as a power law at the phase transition temper¬ 
ature, Tc oc {T — Tpf)~°‘. Consequently, the NMR re¬ 


laxation rate T/^ oc Tc exhibits a sharp divergence at 
T/v. NMR studies of in Ba(Fei_a;Coa;)2As2 and 

BaFe2(Asi_a;Px)2 revealed the presence of spin fluctua¬ 
tions over a broad range of doping and temperature, with 
a quantum phase transition at a critical doping level, Xc, 
that lies close to the maximal 

Several recent experimental studies have reported 
a deviation from the expected power law diver¬ 
gence of as well as stretched exponential be¬ 

havior. In LaFeAsOi-xFa;, Ba(Fei_xRha;)2As2, and 
Ba(Fei_2;COa;)2As2, the characteristic time scale of the 
antiferromagnetic fluctuations grows progressively slower 
over a broad temperature range, the spin-lattice recov¬ 
ery function exhibits stretched exponential behavior, and 
the NMR signal intensity is suppressed (wipeout).^®^^® 
In the case of Ba(Fei_a;Coa;)2As2, also changes 

character in the nematic state, diverging with a crit¬ 
ical exponent S ^ NMR studies at the ®®Co site 
reveal much weaker spin fluctuations near the magnetic 
transition,and ®®Cu site-selective NMR shows a sim¬ 
ilar local suppression of the spin fluctuations on the ®®Cu 
site and neighboring As sites in addition to wipeout the 
NMR signal.^® 

These features point to dynamical inhomogeneity, a 
characteristic of disordered spin glasses indicative of a 
distribution of relaxation rates, in which some fraction of 
the nuclei relax too quickly to be observed.^®’®® Similar 
behavior has been observed in the cluster spin-glass phase 
of the underdoped high Tc cuprates,and charge or¬ 
dering was discovered to be intimately related to the 
®®Cu and ^®®La NQR wipeout in the cuprates.®®’®^ The 
cuprates, however, are doped Mott insulators, and the 
glassy behavior was attributed to intrinsic frustration 
between the competing effects of Coulomb repulsion and 
charge segregation.®®’®® The iron arsenides do not exhibit 
charge ordering and thus a different mechanism must be 
driving the glassy dynamics. In order to investigate the 
glassy behavior in more detail, we have conducted de- 
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X = 0.062 + 0.002 X = 0.066 + 0.007 



FIG. 1. ^®As spectra versus temperature for two different 
doping levels in Ba(Fei_a;Coa;) 2 As 2 measured by sweeping fre¬ 
quency at a constant field of 11.7 T and acquiring echoes for 
the field oriented perpendicular to the c axis. The spectra 
have been normalized to have equal intensities for compari¬ 
son. 



FIG. 2. Magnetization versus recovery time for the ^®As in 
Ba(Fei_j;Goi) 2 As 2 with x = 0.062 at 30K. The solid line is 
the best fit using the protocol described in the text, and the 
dashed and dotted lines are the recovery curves assuming a 
stretched exponential (as described in Ref. 14) or a single 
relaxation time, respectively. 


tailed field, temperature and doping dependent studies 
of both the spin-lattice relaxation rate, , and the 
spin-spin decoherence rate, We extract the tem¬ 

perature dependence of the correlation time, Tc, and find 
that it can be described by Vogel-Fulcher behavior. We 
argue that the dynamical heterogeneity arises because 
the dopants introduce quenched random fields coupling 
to the nematic order. This disorder-induced frustration 
plays a significant role in suppressing antiferromagnetism 
and in the emergence of superconductivity. 


II. RELAXATION MEASUREMENTS 

Single crystals of Ba(Fei_a;Coa,) 2 As 2 and 
Ba(Fei_a;Cua;) 2 As 2 were grown from a FeAs self 
flux and the dopant concentrations were determined 
via wavelength dispersive X-ray spectroscopy (WDS) 
as described in Ref. 27. Multiple WDS measurements 
were made for each batch, and the error bars on the 
concentrations are given by twice the standard deviation 
of these measurements. ^®As (100% abundant, / = 3/2) 
NMR spectra, spin lattice relaxation, and spin echo de¬ 
cays were measured at the central transition (/^ = ±1/2) 
in several different applied fields oriented perpendicular 
to the c axis by acquiring spin echoes using standard 
pulse sequences. Fig. 1 shows representative spectra for 
two different doping levels as a function of temperature. 


A. Relaxation rate distribution 


In order to quantify the distribution of relaxation rates, 
we fit the ^®As magnetization recovery to a distribu¬ 
tion: M{t) = f V{Wi)f{Wit)dWi, where ViyVi) de¬ 
scribes the relaxation rate distribution, and the relax¬ 
ation function f{x) is described below. For a homoge¬ 
neous system V{Wi) is a delta function centered at 
and thus M{t) ^ f{t/Ti). If the distribution has a finite 
width, then the recovery function is more complex, typ¬ 
ically exhibiting stretched behavior. For example, if the 
relaxation function f{x) = e““, then M{t) ^ ^ 

where /3 < 1 is the stretching exponent.^® Previous stud¬ 
ies have reported stretched recovery, however the distri¬ 
bution function for general /3 can only be expressed as 
an infinite series. Here we assume a log-normal distribu¬ 
tion V{Wi) with median and standard devia¬ 

tion (Ti, and fit the magnetization recovery directly. This 
form was chosen because it mimics the distribution for 
a stretched exponential recovery. This approach enables 
us to extract the width of the dynamical distribution of 
the nuclei that contribute to the NMR signal, a quantity 
that sheds important light on the glassy behavior. 

A representative recovery data set with the best fit is 
shown in Fig. 2. The distribution function is given by: 


V{Wi) 


1 I 

- —e ) 

W\a\p2ni: 


( 1 ) 


where ^ and a are variable parameters. We define ^ as 
the median of the distribution, e^, and the standard devi¬ 


ation is given by cti 
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Ba(Feo 938^00 062)2-^^2, Ho — 11.7 T 



Vfi (s’*) 


Ba(Fe,.^Coj2As2 

♦ x = 0.022, Ho = 8.75 T A x = 0.059, Ho = 11.7 T ◄ x = 0.066, Ho = 11.7 T 

• x = 0.058, Ho= 11.7T ■ x = 0.062, Ho = 11.7 T ► x = 0.075, Ho = 3 T 




FIG. 3. (a) Temperature dependence of V{Wi) (normalized by peak height for clarity) and the median ^®As spin-lattice 
relaxation rate, , (markers) for Co-doping with x = 0.062 for Ho || ah. Here is the median of the distribution, 

'P{W\), as described in the text. The probability distribution broadens as temperature is decreased below ~ 100 K. Note the 
bottom axis is a log scale; the skewness of the Log-Normal distribution results in the median falling on the high side of the peak 
(mode), (b) for several Co concentrations as a function of temperature, (c) Standard deviation cti = \/ (Wl) — (VFi)^ of 

the distribution 'P{Wi) for the same samples as a function of temperature in the normal state. Dashed lines in all subfigures 

Q O'/_QA 

indicate structural transition/nematic ordering temperature via bulk measurements reproduced from the literature. ’ 



FIG. 4. Temperature dependence of V{Wi) (normalized by 
peak height for clarity) and the median ^®As spin-lattice re¬ 
laxation rate, Tj”*) (markers) for Co-doping with x = 0.062 
for Ho II ab. The data is identical to that in Fig. 3(a), but is 
plotted on a linear scale. 


function for a spin 3/2 nucleus at the central transition 
is: 


/(x) = Mo 


1 - 2(j) 



( 2 ) 


where x = Wit, t is the time between the initial invert¬ 
ing (or saturating) pulse and the spin echo which samples 


the nuclear magnetization, Mq is the equilibrium nuclear 
magnetization and </> is the inversion fraction. M (t) was 
numerically integrated during fitting using an adaptive 
Gaussian quadrature method and recalculated iteratively 
using a least squares method. The limits for the numeri¬ 
cal integration were chosen to be 10“®s and 10®s, though 
choosing a smaller range when cr is small results in faster 
convergence. This choice of limits was made based on 
the timescale of the NMR experiment. Spins that relax 
faster or slower than this time window will not partici¬ 
pate in the spin echo, and therefore provide natural limits 
of integration. The solid line in Fig. 2 shows the best 
fit determined in this fashion, as well as the best fits as¬ 
suming either a stretched exponential, or a single value 
of using Eq. 2. It is clear that a single uniform re¬ 
laxation rate does not accurately describe the data, but 
both the stretched exponential and the distribution fit 
well. 

Figs. 3(a) and 4 show the temperature dependence 
of the distribution V{Wi), and CTi as a function of 
temperature for Ba(Fei_a;Co 2 ;) 2 As 2 . The data reveal a 
progressive broadening of the distribution below lOOK, as 
well as an increase in both and cri reaching a peak at 
a temperature that coincides with the onset of long-range 
antiferromagnetic order at ■ The peak temperature is 
strongly doping dependent, reflecting the suppression of 
Tjsf with doping concentration. The width ai increases 
by two orders of magnitude, and is also doping depen¬ 
dent. This quantity is a direct measure of the degree of 
dynamical inhomogeneity of the system. Note that at 
low temperatures it is likely that the true width is even 
larger, but we are unable to capture the full distribu¬ 
tion due to signal wipeout. A previous NMR study re- 
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FIG. 5. Field and temperature dependence of the median ) and standard deviation (cti) of the distribution of ^®As 

relaxation rates for Co- and Cu-doped BaFe 2 As 2 . The peak in is strongly field dependent, typical for glassy dynamics. 

ai grows substantially 10^ s“^) below 100 K, reflecting the inhomogeneous relaxation of the nuclei. Dashed lines indicate 
structural transition temperatures via bulk measurements reproduced from the literature. 


vealed that Ba(Fei_a;Coa;) 2 As 2 forms a cluster spin-glass 
state at low temperature below T/v, characterized by a 
distribution of frozen antiferromagnetic domains coexist¬ 
ing with superconductivity.^^ Subsequent neutron scat¬ 
tering work concluded that this cluster spin-glass (or as 
termed by Lu et al. “moment amplitude spin glass”) 
state emerges also in Ba(Fei_a;Nia,) 2 As 2 .^^ The NMR 
data, however indicate that this inhomogeneity begins 
to form at ^ 100 K, well above Tyr, where the spins are 
fluctuating dynamically. This large onset temperature 
suggests that the inhomogeneous fluctuations are unre¬ 
lated to the presence of superconductivity which emerges 
only below Tjv. Furthermore, if the glassy behavior arises 
strictly from disorder and frustration among the spin ex¬ 
change interactions, it is surprising that the inhomogene¬ 
ity would emerge at temperatures well above Tjv, where 
the spin presumably remain uncoupled. 


B. Field and doping dependence 

In order to explore the glassy behavior in more detail, 
we have carried out detailed studies of the field and tem¬ 
perature dependence of V{Wi) as a function of doping in 
both superconducting and non-superconducting samples. 
Changing the magnetic field alters the Larmor frequency, 
enabling one to probe the frequency dependence of the 
slow dynamics. We measured the relaxation in both 
Ba(Fei_xCox) 2 As 2 (up to 30.4 T at the National High 
Magnetic Field Laboratory) and Ba(Fei_a;Cua;) 2 As 2 (up 
to 8.75 T). Fig. 5 shows for several different doping 
concentrations and fields as a function of temperature 
in Ba(Fei_a;Coa;) 2 As 2 and Ba(Fei_a;Cua;) 2 As 2 . Both Co 
and Cu dopants suppress the long range nematic and 
antiferromagnetic ordering, but Cu also suppresses su¬ 
perconductivity to a maximum Ri 2 K, whereas 


reaches a maximum of 23 K in Co-doped samples. 
This enables us to discern whether the glassy behavior is 
connected to the competing superconducting and antifer¬ 
romagnetic ground states.Both systems exhibit qual¬ 
itatively similar glassy behavior, suggesting that its ori¬ 
gin is unrelated to the superconductivity. The maximum 
is suppressed with field, reflecting the fact that 
the relaxation measurement is sampling the fluctuation 
spectrum at a different Larmor frequency. For a hyper- 
fine field h(t), the autocorrelation function is given by 
{h{t)h{0)) = /iQe“‘Ac^ where ho is the root mean square 
value of the field and Tc is the autocorrelation time.^^ In 
this case, the nuclear spin-lattice relaxation rate is: 


w; 




( 3 ) 


where 7 = 7.2919 MHz/T is the ^®As gyromagnetic ra¬ 
tio and ujl = jHa is the NMR Larmor frequency. Note 
that V^Wi) reflects a distribution of both Tc and ho. For 
concreteness, however, we consider only single values of 
these quantities giving rise to the median of the distri¬ 
bution, which is an oversimplification for the real 

system. Eq. 3 shows that reaches a maximum when 
ojlTc = 1 and is equal to T^'^ax = lhl/2Ho. Fig. 6 (a) 
shows T^xnax varies linearly with for various dop¬ 
ings, as expected. The slope of this line gives ho (fit 
values given in Table I), which decreases with dopant 
concentration, and agrees with previous measurements 
in LaFeAsOi-ajF^;.^^ 

Using the measured ho, we proceed to extract Tc. Solv¬ 
ing Eqn. 3 for Tc yields: 


Tc = Wr 




-1 


t; 
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Dopant 

Doping X 

ho (mT) 

offset (s *) 

Co 

0.058 

7.04 ± 0.43 

117.32 ± 24.50 

Co 

0.062 

4.13 ± 0.34 

90.28 ± 11.60 

Cu 

0.040 

4.89 ± 0.07 

48.02 ± 2.77 


TABLE I. Fit parameters extracted for linear fits to 

— l.max / rr—l\ 

1 (-^0 )• 


where the positive sign for the radical arises at low tem¬ 
perature below where Tc ^ and the negative 

sign arises at high temperatures when Tc <C Fig. 

6 (b) presents an Arrhenius plot of TcItco versus T“^, 
where Tco is the high temperature limit of the correlation 
time. The data clearly deviate from linearity, indicating 
that there is not a single activation energy that describes 
the system. The solid black line represents a Vogel- 
Fulcher-Tamman law (tc/tcq = exp{DTK/{T — T^)), 
with D = 0.5(4) and Tk = 25(3)K). This behavior is 
often found in glassy systems, and indicates a ‘fragile’ 
glass, in which the effective activation energy increases 
with decreasing temperature reflecting the collective na¬ 
ture of the fluctuations.^^ Tk represents the temperature 
below which the system becomes trapped in a local min¬ 
imum in free energy at a glass transition temperature. 
In this case, Tk appears to correspond roughly with the 
Neel temperature. However, based on constant field Co¬ 
doping variation fits, this trend appears to break down 
once T/v(NMR) < Tc, where Tk{NMR) is the tempera¬ 
ture at which T^^ reaches a maximum. Below this tem¬ 
perature, the spins are ordered in frozen clusters with 
a broad distribution of sizes.For the Cu-doped sys¬ 
tem, the Tc exhibits more Arrhenius-type behavior. At 
X = 0.04, the peak temperature of Tj“^ is ~ 20 K, which 
agrees with the phase diagram determined via bulk trans¬ 
port and magnetization.^*^ It is unclear why the Cu-doped 
samples differ, but the data suggest that the fluctuations 
are less correlated in this system, which may, in turn, be 
related to the strong suppression of the superconductiv¬ 
ity in this compound. Recent ^'^Cu NMR data suggest 
a strong local effect of the dopants, supporting such an 
interpretation.*^® 

The data in Figs. 5 and 6(b) indicate that V{Wi) is 
slightly modified by the field. In particular, the median 
fluctuation rate and the width cti are suppressed by 
fields up to 30.4 T in the Co-doped sample and 8.75 T 
in the Cu-doped sample. These results suggest that in 
high fields the distribution of domain sizes is narrowing 
and shifting toward smaller domains. Note that because 
of the wipeout effects, these characterizations of the tem¬ 
perature dependence of the glassy behavior may not fully 
capture the behavior of the entire distribution. Since we 
are unable to detect large domains (with correspondingly 
large correlation times Tc) due to wipeout, it is possible 
that the field alters the domain distribution in a manner 
that shifts the weight of the observed distribution to- 



FIG. 6. (a) versus for x = 0.058 and x = 0.062 in 

Ba(Fei_a;Coa,) 2 As 2 and x = 0.04 in Ba(Fei_a;Cua;) 2 As 2 . The 
slope of these data sets reveal the RMS hyperfine field values 
at the As site, and the fit coefficients are detailed in Table I. 
(b) Arrhenius plot of log(Tc/Tco) versus inverse temperature 
for several different fields for Co doping with x = 0.058 and 
Cu doping with x = 0.04. The solid black line shows a Vogel- 
Fulcher-Tamman function, as described in the text. 


wards smaller sizes. Superconductivity in the Co-doped 
samples is also strongly suppressed in these fields, which 
may alter somewhat the domain distribution.^'’ 


C. Spin Echo Decay 

Further evidence for glassy behavior is found in the 
temperature dependence of the '’®As spin-echo decay 
curves. In addition to the increase in Tc/tco and cti, the 
NMR signal intensity gradually becomes suppressed and 
the character of the echo decay changes below 100 K. Fig. 
7(a) shows the echo intensity following a standard echo 
pulse sequence — — for Ba(Fei_a;Coa;) 2 As 2 with 

X = 0.062. The intensity decreases with pulse spacing 
T due to various decoherence effects, including fluctua¬ 
tions of the hyperfine field, h{t), over the course of the 
spin echo experiment. The data have been normalized 
by temperature to account for the Curie susceptibility of 
the nuclei, and clearly reveal the suppression of intensity 
(wipeout) with decreasing temperature.*^ As seen in Fig. 
7(a), the character of the echo decay function crosses over 
from a Gaussian-dominated decay at high temperatures 
to exponential decay below ~ 100 K. This crossover is 
due to the growth of fast spin fluctuations, contributing 
a factor to the echo decay, with W 2 = 

Here is the root mean square of the hyperfine field 
parallel to iJo, in contrast to ho in Eq. 3 which lies 
perpendicular to Hq.^^ Since there is a distribution of 
correlation times Tc as evident from the Tj“* data, we 
fit the echo decay data with the same protocol involving 
a distribution of decoherence rates, IF 2 . The data were 
fit to the function: M(2t) = V{W2)g{2T)dW2, where 

g{T) = / 2 T’ 2 f,g- 2 W 2 T^ Here W 2 is the exponen¬ 

tial component of the spin-spin relaxation rate due to 
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Co-doping X = 0.062, Hq = 8.77 T 



o Co .t = 0.058, Ho = 8.794 Tile A Cox = 0.062, Ho = 3.5 TII afc 
0 Co .t = 0.062, Ho = 8.795 Tile ◄ Co x = 0.062, Ho = 8.77 T II ai 
O CuA- = 0.040,Ho = 4.5Tllab □ Cu a = 0.040, Hj = 6.5 T II ab 
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FIG. 7. Temperature and field dependence of ^®As echo decay for Co- and Cu-doped BaFe 2 As 2 . (a) Echo decay curves 
for X = 0.062, scaled by the nuclear Curie susceptibility. These data are fit (solid lines, see text for details) to extract the 
distribution of spin-spin relaxation rates, 'P{W 2 )- The data were globally fit holding the Gaussian component constant as a 
function of temperature, and by employing the same log-normal distribution form to fit the exponential component, (b) The 
median iT^^) and (c) standard deviation {< 72 ) of the distribution versus temperature for several Co-doped samples, (d) 

(e) (72 versus temperature for the Cu-doped sample, (f) The distribution 'PiW 2 ) and median (markers) versus temperature 
for X = 0.04 at 6.5 T for Ba(Fei_a;Cuj;) 2 As 2 . 


spin-fluctuations, r is the time separating the tt/2 and 
TT pulses of the spin echo sequence, and T 2 G is the tem¬ 
perature independent Gaussian component of the spin- 
spin relaxation. At high temperatures the echo decay 
has a Gaussian form, which reflects the complex direct 
and indirect couplings between the like As nuclei. We do 
not expect this component to change with temperature, 
whereas the growth of spin fluctuations at low temper¬ 
ature will affect Each temperature dependent 

data set was fit globally with a temperature-independent 
T 2 G to achieve the best fit to all temperatures. This 
global analysis was confirmed by individually fitting the 
data set at each temperature, results of which show no 
trend in T 2 G as a function of temperature. 

The data in Fig. 7(f) shows 7^(W2), panels (b) and 
(d) show the median and panels (c) and (e) show 
the standard deviation 172 for several doping levels and 
dopants as a function of temperature. The temperature 
dependence of agrees qualitatively with the corre¬ 
lation times extracted from the data seen in Fig. 

6(b). T 2 ^ = therefore we expect a monotonic 

increase of with decreasing temperature. Surpris¬ 
ingly, the width 02 of this distribution differs from ui 
extracted from the spin-lattice relaxation data, and ex¬ 
hibits a downturn below Tm- Note, however, that 'P{W 2 ) 
is cut off at large W 2 by the finite detection window of the 
NMR experiment, which is the primary cause of signal 
wipeout.As a result, the measured width a 2 is reduced 
as the majority of the distribution shifts outside of the 


detection window at low temperature. 


III. DISCUSSION 
A. Missing signal 

It is clear from Fig. 7 that signal wipeout of up to 80% 
is present, consistent with previous measurements of the 
spectral intensity in these samples, which raises the ques¬ 
tion of where the missing signal has gone.^^ The system 
is either dynamically or spectrally inhomogeneous. In 
our experiments we find no significant broadening of the 
spectra in the paramagnetic state, as shown in Fig. 1. It 
is possible that the distribution is such that a large frac¬ 
tion of the nuclei resonate outside of this window, but the 
internal field in the ordered state is small in this range of 
dopings and the spectral shift for this field orientation is 
minimal.^® The spin lattice relaxation was measured at 
the peak of this resonance, and it is possible that not all 
of the nuclei were inverted by the radiofrequency pulses. 
It is more likely, however, that the missing signal arises 
from dynamical heterogeneity, given the broad distribu¬ 
tion of relaxation rates that we observe. The missing 
signal in this case arises from nuclei that are located in 
an environment with a sufficiently long Tc such that they 
recover to equilibrium before they can contribute to the 
spin-echo signal. It is important to note that the distribu¬ 
tions shown in Fig. 3, 4 and 7 are representative only of 
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the nuclei that are actually contributing to the signal. In 
fact the true distributions are likely to be much broader 
than what we are able to measure, as a significant portion 
of the nuclei experience even faster relaxation rates. 


B. Glassy nematic fluctuations 

The inhomogeneous fluctuation distribution cannot be 
understood simply in terms of critical slowing down of 
the spin degrees of freedom. The spin fluctuations are 
not averaged out spatially, implying the existence of mul¬ 
tiple local domains of characteristic size ^ ^ Tc- Figure 
8 (a) summarizes the doping dependence of the width, 
ai{x,T), of the inhomogeneous distribution, where cti is 
related to the distribution of domain sizes. For the parent 
compound BaFe 2 As 2 we find that the recovery fits best 
to a single component of relaxation for all temperatures, 
so the system is dynamically homogeneous, cti remains 
small for the lightly Co-doped regime; however near opti¬ 
mal doping it becomes a strong function of temperature, 
reflecting a large dynamical inhomogeneity both in the 
Co and Cu-doped crystals. 

A likely origin for this inhomogeneity is nematic fluc¬ 
tuations associated with the proximate tetragonal-to- 
orthorhombic structural phase transition. The doping 
and temperature trends exhibited by <7i{x,T) shown in 
Fig. 8 (a) closely resemble the behavior of the static ne¬ 
matic susceptibility, Chu et al. have found that 

Xn{T) exhibits Curie-Weiss behavior, with Weiss temper¬ 
ature 6 that vanishes at the critical doping of Xc = 0.07 
for the Co-doped system. The fluctuation-dissipation 
theorem implies that because of the large susceptibility, 
there are also significant thermal fluctuations of the ne¬ 
matic order. In other words, even though there is no long- 
range nematic order, local orthorhombic distortions con¬ 
tinue to fluctuate well above the ordering temperature. 
Because the spins are strongly coupled to the nematicity, 
these nematic fluctuations will drive spin fluctuations, 
which in turn couple to the nuclei via the hyperfine inter¬ 
action to influence nuclear spin-lattice relaxation. In fact, 
scales with shear modulus in this phase, reflecting 
the fact that both quantities are probing the dynamics 
of the nematic fluctuations.^^ 

The glassy inhomogeneous nature of the fluctuations, 
therefore, probably reflects a property of the nematic 
fluctuations. Because the nematic order has Ising sym¬ 
metry and breaks spatial symmetry, it is highly sensitive 
to quenched random impurities and is prone to exhibit 
glassy behavior.^® The theory of electronic nematic or¬ 
der and the role of disorder is well established in the 
context of the cuprates,and more recently in the 
context of the iron pnictides.^®’^® The dopant atoms may 
provide a random field potential for nematic order that 
suppresses the phase transition temperature and gives 
rise to a distribution of frustrated nematic domains with 
different fluctuation rates, as illustrated in Fig. 8 (b). 
With increasing dopant concentration, the nematic or- 



X in Ba(Fe|,^-CUj-) 2 As 2 x in Ba(Fe|_,jCo,.) 2 As 2 



FIG. 8. (a) Temperature vs. doping phase diagram of Cu- 
and Co-doped BaFe 2 As 2 . Markers have been reproduced from 
bulk measurements in the literature, and solid lines are a 
guide to the The color scale overlay shows the 

standard deviation ui = \/~ (IFi)^ for the distribu¬ 
tion 'P{Wi), characterizing the degree of inhomogeneity of 
the NMR spin-lattice relaxation rate, (b) Schematic of local 
nematic domains, indicating directions of Fe spin (arrows). 
The tetragonal and orthorhombic unit cell axes are shown. 
The local nematicity is oriented along the ellipses. 


dering transition is gradually suppressed. In the disor¬ 
dered phase, there are fluctuating patches in which C 4 
symmetry is temporarily and locally broken, but there is 
no long range or static order. These fluctuating patches, 
however, exhibit a broad range of sizes and fluctuation 
times. The inhomogeneity we observe reflects the distri¬ 
bution of these patches. The NMR data indicate that the 
nematic fluctuations and distribution of domains persist 
up to ~ 100 K, as shown in Fig. 8 (a). The local autocor¬ 
relation time of the domains, Tc, is proportional to the 
domain size, thus the width of the distribution of domain 
sizes grows up to two orders of magnitude by the onset 
of long range nematic order. This scenario provides a 
natural explanation for the large Xn as well as the broad 
distribution of relaxation times observed in our NMR ex¬ 
periments. Further, it explains the similarity of the phase 
diagram of both electron and hole-doped systems, as well 
as the isovalent BaFe 2 (Asi_ 2 ;Pa ;)2 system. 

The temperature-pressure phase diagram of the stoi¬ 
chiometric parent compound also exhibits a suppression 
of antiferromagnetism and emergence of superconductiv¬ 
ity without the presence of dopants;^^’^® however, natural 
lattice defects may provide a source of quenched disor- 





















der that could be amplified by non-hydrostatic pressure. 
On the other hand, non-isovalent dopants clearly play 
a role in tuning the density of states, as revealed by a 
recent study of simultaneous hole- and electron-doping 
in Baia;K 3 ;Fei. 86 Coo.i 4 As 2 demonstrating that that the 
magnetic state can be partially recovered by compensat¬ 
ing the carrier concentration.^® Thus both disorder and 
tuning the density of states appear to be important pa¬ 
rameters controlling the phase diagram. 

IV. CONCLUSION 

In summary, the glassy behavior we observe in the dy¬ 
namics reveal a highly inhomogeneous system in a re¬ 
gion of the phase diagram that is nominally a homo¬ 
geneous disordered paramagnetic phase. The NMR re¬ 
sponse probes the Fe spins through the hyperfine cou¬ 
pling, but it is the nematicity that drives the response of 
the system. The disorder introduced by the dopants gen¬ 
erates random strain fields, which couple to the nematic¬ 
ity and may contribute to the suppression of the nematic 
ordering temperature. The nematic order parameter de¬ 
velops a complex fluctuating spatial landscape, with var¬ 
ious domain sizes. Future measurements under uniaxial 
strain may significantly suppress the width of the distri¬ 
bution, and will provide an important avenue to investi¬ 


gate the dynamics in the glassy phase. NMR studies of 
the dynamics under pressures up to 10-15 GPa in stoi¬ 
chiometric samples will also help to elucidate the role of 
disorder in suppressing the nematic phase. 
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